Nature of the insulating phases in the half-filled ionic Hubbard model 
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We investigate the ground-state phase diagram of the one-dimensional "ionic" Hubbard model with 
an alternating periodic potential at half-filling by numerical diagonalization of finite systems with 
the Lanczos and density matrix renormalization group (DMRG) methods. We identify an insulator- 
insulator phase transition from a band to a correlated insulator with simultaneous charge and bond- 
charge order. The transition point is characterized by the vanishing of the optical excitation gap 
while simultaneously the charge and spin gaps remain finite and equal. Indications for a possible 
second transition into a Mott-insulator phase are discussed. 



PACS numbers: 71.10.-w, 71.10.Fd, 71.10.Hf, 71.27.-Ha, 71.30.-hh 



I. INTRODUCTION 

For more than two decades the correlation induced 
metal-insulator transition and its characteristics has been 
one of the challenging problems in condensed matter 
physics [1]. This metal-insulator transition is often ac- 
companied by a symmetry breaking and the develop- 
ment of long range order [2]. In one dimension this or- 
dering can only be related to the breaking of a discrete 
symmetry. Examples include commensurate charge den- 
sity waves (CDWs) and Peierls dimerization (bond-order 
wave, BOW) phenomena. In contrast, the transition into 
the Mott insulating (MI) phase in one dimension is not 
connected with the breaking of a discrete symmetry [3]. 
In the MI phase the gapped charge degrees of freedom 
are uniformly distributed in the system, while the gap- 
less spin degrees of freedom are described by an effective 
S = 1/2 Heisenberg chain [4]. 

Due to the different symmetries of the CDW, BOW, 
and MI phases it is natural to expect that these phases 
are mutually exclusive. The extended Hubbard model at 
half-filling with an on-site {U) and a nearest neighbor {V) 
Coulomb repulsion provides a prominent example with a 
transition from a MI to a CDW insulator in the vicinity 
of the U = 2V line in the phase diagram [5] . Remarkably, 
the transition at weak coupling may involve an intermedi- 
ate BOW phase [6,7]. A similar phase-diagram structure 
was recently also discovered for the Holstein Hubbard 
model [8]. The tendency towards BOW order is even 
more profound for the Hubbard model with an explicit 
bond-charge coupling where the CDW and MI phases are 
often separated by a long range ordered Peierls dimerized 
phase [9]. 

In recent years particular attention has been given to 
another example for an extension of the Hubbard model 
which includes a staggered potential term [10-18]. The 
corresponding Hamiltonian has been named the "ionic 
Hubbard model" (IHM); in one dimension it is given by 



(1) 



where creates an electron on site i with spin a and 

"io- ~ '^la^ia- ^ potential energy difference be- 

tween neighboring sites, and 5 a Peierls modulation of 
the hopping amplitude t. In the limit A = 5 = 0, Eq. 
(1) reduces to the ordinary Hubbard model, the limit 
A = and (5 > is called the Peierls-Hubbard model, 
and the limit A > and 5 = is usually referred to as 
the IHM. In the following, we will focus mainly on the ef- 
fect of the on-site modulation A, so we implicitly assume 
(5 = except where stated otherwise. 

The IHM was first proposed and discussed almost 20 
years ago in the context of organic mixed-stack charge- 
transfer crystals with alternating donor {D) and accep- 
tor {A) molecules (. . . D+pA'P . . .) [19,10]. These stacks 
form quasi-lD insulating chains, and at room tempera- 
ture and ambient pressure are either mostly ionic (p w 1) 
or mostly neutral {p « 0) [19]. However, several sys- 
tems undergo a reversible neutral to ionic phase tran- 
sition i.e. a discontinuous jump in the ionicity p upon 
changing temperature or pressure [20]. Later the IHM 
has been used in a similar context to describe the ferro- 
electric transition in perovskite materials such as BaTiOa 
[21,22] or KNbOg [23]. 

The very presence of at least on transition in the 
ground state phase diagram of the half-filled IHM model 
is easily traced by starting from the atomic limit [24,13]. 
For t = 0, it is obvious that at [/ < A the ground state 
of the IHM has two electrons on the odd sites, and no 
electrons on the even sites corresponding to CDW or- 
der with maximum amplitude. On the other hand, for 
t/ > A each site is occupied by one electron with infinite 
spin degeneracy. Thus, for i = a transition occurs at 
a critical value Uc = A. This transition is expected to 
persist for finite hopping amplitudes t > 0. 

A renewal of interest in the IHM started with the 
bosonization analysis of Fabrizio et al. (FGN) [12], where 
a two-transition scenario for the ground-state phase dia- 
gram of the ID IHM was proposed. The key features of 
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the FGN theory are the presence of an Ising-type tran- 
sition from a CDW band-insulator phase at [/ < L/^^ 
into a BOW phase at U^/^ < U < U^p, and a continu- 
ous Kosterlitz- Thouless like transition into a MI phase at 
U > Ugp. In this scenario the charge gap vanishes only 
at [/ = U^f^ and the system might be "metallic" at this 
point. The second transition at C/ = C/jp is connected 
with the closing of the spin gap, which is finite for all 
U < and vanishes for U > Ug^. Thus, the bosoniza- 
tion phase diagram essentially supports the "exclusion 
principle" of the ground states. 

Later on various attempts based on numerical tools 
have been performed to verify the FGN phase diagram 
for the half-filled IHM. In particular, exact diagonal- 
ization [13,18], valence bond techniques [15], quantum 
Monte Carlo [14], and DMRG [16,17] was used. Unfor- 
tunately, conflicting results have so far been reported in 
these studies regarding the nature of the transition and 
the insulating phases, the possibility of two rather than 
one critical point, or the appearance of BOW order. 

Given the numerous unresolved issues we reinvestigate 
the ground-state properties of the IHM using the ex- 
act diagonalization Lanczos technique and the DMRG 
method. We verify the presence of at least one transition 
at a critical coupling C/c(A) from a band-insulator (BI) 
to a correlated insulator (CI) phase. On finite systems 
the transition originates from a ground-state level cross- 
ing with a change of the site-parity eigenvalue, which 
implies the vanishing of the optical excitation gap at Uc] 
Our DMRG results show that the spin and charge gaps 
remain nevertheless finite and equal at the transition. 
Above Uc the charge and spin gaps split, the charge gap 
increases, while the spin gap decreases and we identify 
long-range BOW order with a spontaneous site-inversion 
symmetry breaking. The existence of a second transition 
is not unambiguously resolved within the accuracy of our 
DMRG data. Yet, the scaling of the BOW order parame- 
ter changes qualitatively with increasing U indicative for 
a possible second smooth transition point where the spin 
gap closes. 

We show that aXU < Ud^) the CDW-band insulator 
phase is realized. In this phase BOW and spin density 
wave (SDW) correlations are strongly suppressed, and 
the spin and charge gaps are equal and finite. The charac- 
teristic feature of the ground-state phases for U > ?7c(A) 
is the coexistence of long range CDW order with either 
a long range BOW or algebraically decaying BOW and 
SDW correlations. 



II. SYMMETRY ANALYSIS 

A good starting point for understanding the existence 
of a phase transition in the IHM is to study the symmetry 
of the model manifestly seen in the limiting cases U <C 
A,t and [/ 3> A, t. The IHM is invariant with respect 
to inversion at a site and translation by two lattice sites. 



If we denote the site inversion operator by P, defined 
through 

p4p+=cL,. for i = 0, 1 , (2) 

and Tj for a translation by j sites, then any nondegener- 
ate eigenstate jf/'n) of H must obey P\tpn) = ±|V'n) and 
T2\'ipn) = IV'n). Because [H,Ti] ^ 0, any non-degenerate 
eigenstate jf/'n) of H is not an eigenstate of Ti. 

For the half-filled Hubbard model {A = S = 0) the 
ground state has P = +1 only for U = 0, and P = — 1 for 
any U > [13]. However, in the IHM the phase transition 
from a renormalized BI to a CI occurs at some finite 
Z7c > 0. This suggests that the parity of the ground state 
remains even not only for U = 0, but for all U < Uc- At 
Uc, a ground-state level crossing occurs on finite chains, 
as confirmed by exact diagonalization studies (see below), 
connected with a site-parity change. 

For U = the ground state at half-fiUing is a CDW-BI. 
The alternating potential defines two sublattices, dou- 
bling the unit cell and opening up a band gap A for 
?7 = at fc = ±7r/2. The elementary spectrum con- 
sists of particle-hole excitations over the band gap. The 
charge (Ac) and spin {As) excitation gaps are equal: 
Ac = As = A. We consider a system to be a BI when 
the criterion A5 = Ac holds, where the spin and the 
charge gaps arc given by 

As = EoiN = L,S, = 1) 
-Eo{N = L,S,=0), 

Ac = Eo{N = L + l,S, = l/2) 
+ Eo{N = L-l,S, = l/2) 
-2Eo{N = L,S, = 0), (3) 

respectively. Eo{N, Sz) is the ground-state energy, L the 
system length, N the number of electrons, and Sz the 
^-component of the total spin. As we show below the 
BI phase is realized in the ground state of the IHM at 
U<Uc. 

In the strong-coupling limit U A, the low-energy 
physics of the IHM is described by the following effective 
Heisenberg spin model [10] 
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In Eq. (4) the exchange couplings are given by 
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where x = A/U. This result (4) implies that in the 
strong-coupling limit of the IHM the low-energy physics 
is qualitatively similar to that of the Hubbard model, 
with modified exchange coupling constants J and J'. For 
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next-nearest neighbor couplings J' < 0.24J the spin gap 
vanishes [25]. The couphng constants (5) satisfy the con- 
dition J' < 0.24J at least for U > 3.6t for A < i and 
U > 3.6A for A > i. 

The effective spin model (4) is invariant with respect to 
translations by one lattice spacing, whereas the original 
IHM is invariant only with respect to translations by two 
lattice spacings. However, the doubling of the unit cell 
is ensured due to the charge degrees of freedom. For the 
standard Hubbard model at arbitrary U the number 
of doubly occupied sites D in the ground state is finite. 
The exact Bcthe-ansatz solution tells that D scales as 
{t/UY in the strong coupling limit and is given by [26] 

D = ^ NA{t/Uni + O {{t/Uf)], (6) 

i 

where A = 4 In 2. Contrary to the Hubbard model, where 
doublons are equally distributed on all sites of the system, 
the non-equivalence of sites in the IHM leads to different 
probabilities for finding a doublon on even or odd sites. 
Since doublons are spin singlets, their distribution is not 
influenced by the spin fluctuations. Since the energy of 
doublons on even and odd sites differ in A and assuming 
the scaling for the density of doublons as in Eq. (6) , one 
easily obtains for the amplitude of the ionicity induced 
CDW in the strong coupling limit 
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C/2 (1 - x2)2 



[l + 0{it/ur)] (7) 



where Ai is a constant of order unity. Thus, although 
the effective spin Hamiltonian has a higher symmetry 
than the original model from which it was derived, the 
translational symmetry of the IHM is recovered due to 
the long range CDW pattern arising from the staggered 
doublon and holon distribution. 



III. EXACT DIAGONALIZATION RESULTS 

In order to explore the nature of the spectrum and 
the phase transition, wc have diagonalized numerically 
small systems by the Lanczos method [27] similarly to 
earlier exact diagonalization calculations [11,13]. The 
energies of the few lowest eigcnstatcs were obtained for 
finite chains with L = 4n and periodic boundary condi- 
tions or L = 4n + 2 with antiperiodic boundary condi- 
tions. 

We first analyze short chains: for chain lengths L < 16 
finitc;-sizc effects do not change; the qualitative behavior 
discussed below. In Fig. 1, the lowest eigenenergies of 
the IHM for A = 0.5t, L = 8 and periodic boundary 
condition are shown as a function of U. At U = 1.3t, a 
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FIG. 1. Lowest energy eigenvalues of the IHM at half-filling 
for L = 8 sites, periodic boundary conditions and A = 0.5t. 



level crossing of the two lowest eigenstates occurs. A non- 
degenerate eigenstate of the IHM has a well defined site- 
parity, so a ground-state level-crossing transition neces- 
sarily corresponds to a change of the site-parity eigen- 
value. 

For J7 = 0, the IHM is easily diagonalized in momen- 
tum space by introducing fermionic creation operators 
Tfcab with a band index 6 = 1,2 for the lower and up- 
per bands, respectively, with the dispersion Ei/2{k) = 
±A/4cos2(fc) + (A/4)2 for momenta -tt/2 < k < tt/2. 
For [7 = the first two degenerate excited states at 
half-filling always have negative site parity, because the 
ground state has P = +1, and the operator 7g<j27g(7i with 
q = 7r/2 obeys 



-^^70(7270(71 ~ 7g<727(7(7l-P 



(8) 



The first two excited states shown in Fig. 1 are the spin 
singlet {S = 0, Sz = 0) and triplet excitations {S = 1, 
Sz =0), created from the ground state by applying the 
operators 
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(7gT27,|l + 7g427g4i) , 



(9) 



respectively. Thus both excited states have total momen- 
tum ktot = and negative site parity. For U > 0, these 
degenerate excited states split in energy. Exact diagonal- 
ization of finite IHM rings therefore identifies one critical 
Uc > 0, separating a BI with P = +1 at U < Uc from a 
CIwithP = -lfor(7>C/c. 
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IV. DMRG RESULTS 

A. Excitation gaps 

In order to access the transition scenario in the long 
chain-length limit, we have studied chains up to L=512 
using the DMRG method [28-30]. The fact that the 
transition at Uc is connected to a change in inversion 
symmetry requires some caution when open boundary 
conditions (OBC) are used in DMRG studies. For OBC 
and L = 2n the IHM is not reflection symmetric at any 
site. Thus, the ground state does not have a well defined 
site parity, and the level-crossing transition is absent. To 
overcome this problem, one might try to use chains with 
OBC and an odd number of sites L = 2rt + 1, since the 
Hamiltonian in this case is reflection symmetric with re- 
spect to the site ic in the center of the chain, and a site 
inversion operator is well defined by 

p4..P^ = ^ = 0,...,L-l. (10) 

To test whether this is an improved choice we have cal- 
culated the site parity of the ground state for [/ = 
analytically for different chain lengths L = 2n + \ and 
found 

PlV-o) = (-l)"|V'o) . (11) 

On the other hand, if one extends the idea of Gidopoulos 
et al. [13] for the determination of the site parity to chains 
with L = 2n + 1 for U ^ t, one obtains 

p|^o) = (-i)El:i"] i^o) = (-i)"iv^o) . (12) 

Thus, the parity eigenvalue of the ground state is the 
same at [/ = and U ^ t for a given chain length, 
and no level crossing occurs. Due to the fact that the 
sharp transition at a well defined Uc does not exist in the 
finite-chain results for OBC, the extrapolation is a rather 
subtle problem, since a sharp transition feature has to be 
identified from the extrapolation of smooth curves. This 
requires the use of quite long chains in the critical region. 

In Fig. 2 extrapolated results are shown for the spin 
and charge gaps As and Ac, respectively. Calcula- 
tions were performed with OBC for chains of lengths 
L = {30,40,50,60}, and additionally up to L = 512 
in the transition region around the estimated Uc- We 
assume a scaling behavior of Ac and A5 of the form [32] 

A,(L) = Ar + ^ + §, (13) 

where i € {S,C}. The extrapolation for L — s- cx) is 
then performed by fitting this polynomial in l/L to the 
calculated finite-chain results. We note that different 
finite-size scaling formulas were proposed in the litera- 
ture mainly when periodic or antiperiodic boundary con- 
ditions were used [13]. 
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FIG. 2. Results for the spin (As) and charge (Ac) gaps 
of the IHM at half-filling with A = O.bt as a function of 
U. Energies were obtained by DMRG calculations on open 
chains with L = {30, 40, 50, 60} (main plot) and up to i = 512 
(inset), and extrapolated to the limit of infinite chain length. 

As can be seen from the main plot in Fig. 2, extrapo- 
lating the results for L — {30, 40, 50, 60} does indeed not 
give a sharp transition behavior. As illustrated in the 
inset, adding results for L up to 512 in the critical region 
changes the picture considerably. Within numerical ac- 
curacy the charge and spin gaps remain equal and finite 
up to a critical Uc ~ 2. It, where a sharp kink for Ac is 
observed. Importantly, Ac does not close at the critical 
point. We emphasize that the magnitude of Ac and A5 
at the transition point is sufficiently larger than our nu- 
merical uncertainty in the finite size scaling analysis and 
therefore allows for a safe conclusion. Ac — A5 > at 
the transition is in fact not in confiict with an underly- 
ing ground-state level crossing. If the ground states of 
the different site-parity sectors become degenerate, the 
only rigorous consequence is the closing of the optical 
excitation gap. The selection rules for optical excitations 
allow only for transitions between states of different site- 
parity. Furthermore, optical transitions occur within the 
same particle number sector. The optical gap is therefore 
by definition distinct from the charge gap Eq. (2) which 
involves the removal or the addition of a particle. The 
critical point Uc of the IHM has the remarkable peculiar- 
ity that the optical gap closes while Ac remains finite. 
Above Uc the charge and spin gaps split indicating that 
the corresponding insulating phase is no longer a BI. A5 
continuously decreases with increasing U and becomes 
unresolvably small above U ~ 2.5t within the achievable 
numerical accuracy. 

The result, that Ac and As remain finite and equal at 
the transition is in agreement with the data obtained by 
Qin et al. [17]. These authors performed DMRG calcu- 
lations for the IHM with A — 0.6t for open chains up to 
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L = 600 sites. They observed a surprising non-monotonic 
scaling behavior of A5 with L for values of U close to the 
critical Uc, i-e. for chain lengths L > 300 A5 started to 
increase again. It remains unclear whether this is due 
to loss of DMRG accuracy with increasing chain lengths 
and keeping a fixed number of states in the DMRG algo- 
rithm. In contrast, our data always show a monotonous 
scaling with 1/L. DMRG calculations for the IHM with 
A = t have also been performed by Takada and Kido 
for chains up to L = 400 sites [16]. The authors in- 
terpret their results in the region close to Uc in favour 
of a two-transition scenario similar to that of FGN [12]. 
However, as we will show below the accuracy of the cur- 
rently available DMRG data is not enough to provide a 
stringent argument in favour of this interpretation. 

For comparison we show in Fig. 3 the spin and charge 
gaps versus U in the Peierls Hubbard model. As we ob- 
serve this model is distinctly different from the IHM. This 
is also a band-insulator at U — 0, but in contrast to 
the IHM has > As > for any value U > 0, i.e. 
the phase transition from the Peierls band-insulator to 
the correlated insulator occurs at Uc = 0- So although 
the Peierls and the ionic BI for U — similarly possess 
an excitation gap at the Brillouin zone boundary, apply- 
ing a Coulomb U leads to distinctly different behavior in 
both cases. The origin of the different behaviors must 
be traced to the fact that the Hubbard interaction and 
the ionic potential compete locally on each site, while 
the Peierls modulation of the hopping amplitude tends to 
move charge to the bonds between sites, thereby avoiding 
conflict with the Hubbard term. 




FIG. 3. Results for the spin (Ag) and charge (Ac) gaps 
vs. [/ of the Peierls-Hubbard model at half-filling with a 
modulation of the hopping amplitude 5 = 0.5. Energies 
were obtained by DMRG calculations on open chains with 
L = {30,40, 50,60}. 

The spin-Peierls physics of the Peierls Hubbard model 
at large U evolves smoothly with decreasing U into the 
physics of a spin-gapped CI-BOW state in the weak- 



coupling limit, which is characterized by long-range stag- 
gered bond-density correlations 

i 

m = E + H.c) . (15) 



For U ^ t the low-energy physics of the Peierls Hub- 
bard model is described by the spin-Peierls Heisenberg 
Hamiltonian with a staggered exchange interaction and 
a dimerization induced spin gap [31]. 



B. Correlation functions 

The important question remains about the nature of 
the insulating phase of the IHM for U > Uc- To further 
analyze the BI and CI phases below and above Uc, we 
have evaluated site- and bond-charge distribution func- 
tions as well as spin-spin correlation functions. In 




FIG. 4. Electron density distribution in the ground state 
of the IHM for A = 0.5t and U = 0.8t (diamonds) and U = it 
(circles). Results were obtained by DMRG calculations on an 
open L = 32 chain. 

Fig. 4 we show the charge distribution (0|(n(r) — 1)|0) 
{n{r) ~ n.r.^+nr.i) in the ground state |0) for the IHM at 
U = 0.8t < Uc and [/ = 4t > C/c for a L = 32 chain. The 
alternating pattern in the density distribution is well pro- 
nounced not only in the BI but also in the CI phase far be- 
yond the critical point at U ^ Uc- For the L = 32 chain 
the CDW is well established at distances I ^ L/2 even at 
U — At and its amplitude remains almost unchanged in 
the finite-size scaling analysis (see upper inset in Fig. 5). 
The main plot in Fig. 5 shows that the staggered com- 
ponent of the charge density decreases smoothly with in- 
creasing U - Close above the transition point near U — 2t 
one observes an anomaly, i.e. a slight decrease of the 
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FIG. 5. Staggered charge density component vs. U, for 
A — 0.5t and L — 512 (main plot). Its scaling behavior, for 
U — (diamonds), U = 2.1t (circles) and U — 4.0t (triangles) 
is shown in the inset. 

CDW amplitude accompanied by a change in curvature; 
this anomaly can be clearly identified in the enlargement 
shown in the lower inset. The anomaly will find a natural 
explanation in the discussion below. Our numerical data 
show that the alternating pattern in the electron density 
distribution in the IHM remains for arbitrary finite U. 
Thus, the ionic potential enforces long range CDW order 
for all interaction strengths. 

Fig. 6 shows the DMRG results for the spin-spin cor- 
relation function (0|5^(L/2)S'^(L/2 + r)|0) for the IHM 
a.t U = 0.8t < Uc and U = At > Uc in. comparison with 
the Hubbard chain. In the BI phase at U = 0.8t the 
SDW correlations are quickly suppressed after a few lat- 
tice spacings. At U — At the amplitude of the SDW 
correlations in the CI phase of the IHM is slightly re- 
duced in comparison to the Hubbard model at the same 
value of U. However, the large distance behaviors of the 
spin correlations in the CI phase of the IHM and the MI 
phase of the Hubbard model are quite similar and be- 
come almost indistinguishable (see the data in the inset 
for U = At). On the other hand, the long distance be- 
havior of the spin-spin correlation function at C/ = 2.1t 
and A = 0.5t, i.e. close above the transition, manifestly 
supports the finiteness of the spin gap (see inset) . This 
may be viewed as an indication for existance of two dif- 
ferent phases above Uc- But we cautiously point out that 
it is hard to judge on the persistence or vanishing of As 
far above Uc in the CI phase from the finite chain spin 
correlators alone. 

To address the BOW ordering tendencies in the CI 
phase we have calculated the ground-state distribution 
of the bond-charge density Eq. (15). Fig. 7 shows the 
results of the DMRG calculations for the L ^ 32 IHM 
and the Hubbard chain at U = 0.8t and U = 2.6t. The 
boundary effect for an open chain is strong and leads to 
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FIG. 6. Spin-spin correlation function in the ground state 
of the IHM for A = Q.5t (open symbols) and the Hubbard 
model (A = 0) (full symbols) at U — 0.8t (diamonds) and 
U = At (circles). Chain length L — 32 (main plot). Inset: 
Scahng of the spin-spin correlation function at r = L/4 -I- 1 
for U — 2.1t (up and down triangles) and U = At (circles) 
for the IHM (open symbols) and the Hubbard model (full 
symbols) 

a modulation of the bond density already for the pure 
Hubbard model. Interestingly, the same behavior was 
also observed previously for the bond expectation value 
{Si ■ Si+i) in open antiferromagnetic Heisenberg spin- 
1/2 chains [29]. The detailed comparison with the Hub- 
bard chain in Fig. 7 reveals, that the ionic potential 
leads to a reduction of the bond-density oscillations at 
U < Uc, while in the CI phase at U > Uc their amplitude 
slightly increases. The enhancement of BOW correla- 
tions above Uc must simultaneously weaken the CDW 
amplitude. This naturally explains the slight downward 
curvature near Uc in the staggered charge density com- 
ponent shown in Fig. 5. 

In order to explore the possibility towards true long 
range BOW ordering in the IHM above Uc we plot in 
Fig. 8 the staggered bond-density versus U in the center 
of long, open chains with L ~ 256 and L = 512. In 
the BI phase this quantity is essentially zero. At the 
transition point the staggered component of the bond 
density increases rapidly, and on further increasing U 
it starts to decrease smoothly. However, the staggered 
bond density remains finite for any U > Uc on these long 
but finite chains and vanishes only in the limit [/ — > oo. 

Naturally it is necessary to perform a finite size scal- 
ing analysis for this quantity. Fig. 9 shows its chain 
length scaling behavior for [/ = 0, t/ = 2. It, and U — At. 
Two conclusions can be drawn from these results: In 
the absence of the interaction the staggered bond density 
clearly extrapolates to zero - as expected for a conven- 
tional BI. In the CI phase close above the transition the 
upward curvature of the staggered bond density vs. 1/L 
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FIG. 7. Bond-charge density of the IHM for A = 0.5t (open 
symbols) and the Hubbard model (full symbols) at U = 0.8t 
(diamonds) and U = 2.6t (circles). 

points to a finite value in the infinite chain length limit, 
i.e. long range BOW order. However, the scaling behav- 
ior changes in the CI phase far above the transition point. 
Interestingly, the scaling behavior in the IHM far above 
Uc starts to reseamble the results for the pure Hubbard 
model, for which the staggered bond density has to vanish 
in the thermodynamic limit. The qualitative change in 
the scaling behavior of the staggered bond density may 
again be viewed indicative for a possible second phase 
transition. Summarizing the results for the SDW and 
BOW correlations we conclude that there is evidence for 
two phases for U > Uc- Close above Uc long range BOW 
order develops in the ground state of IHM while far above 
the transition point the coorrelation functions become al- 
most identical to those of the Hubbard model. Yet, long 
range CDW order exists for all U. A precise location of a 
second transition point is however not possible from the 
currently awailablc DMRG data. 

V. CONCLUSIONS 

From the finite chain DMRG studies and finite size 
scaling analyses we draw the following conclusions for the 
ground-state phase diagram of the IHM: the ionic poten- 
tial leads to long range CDW order for all interaction 
strengths. The data resolve one transition point from 
the BI to the CI phase. Remarkably, at the transition 
Ac — A5 and both remain finite. Close above the tran- 
sition, i.e. for < {U/Uc) — 1 <C 1, we identify a clear sig- 
nal for long range BOW order. With increasing U above 
Uc the finite size scaling behavior of the staggered bond 
density and spin-spin correlation function changes quali- 
tatively and approaches the scaling behavior of the Hub- 
bard model. With the current chain length and DMRG 
accuracy limitations it is not possible to precisely iden- 




u/t 

FIG. 8. Staggered bond-density component vs. U near the 
center of the L = 256 (circle, dotted line) and L = 512 (dia- 
monds, solid line) chain. 



tify and locate the second transition point. The phase 
with BOW order necessarily has a finite spin excitation 
gap. If BOW seizes to exist above a second critical value 
of U, the spin gap has to vanish simultaneously leading 
to identical large distance decays of SDW and BOW cor- 
relations. 

The insulator-insulator transition at Uc{A) on finite 
periodic chains results from a ground-state level crossing 
of the two site-parity sectors. The optical excitation gap 
therefore has to vanish at Uc, remarkably, the DMRG 
data reveal that at the critical point Ac ans A5 remain 
both finite and equal. This itself clearly indicates the ex- 
istanse of a CI phase with Ac > A5 > which originates 
from the appearance of long range staggered bond den- 
sity order. The distinction between the optical and the 
charge gap is therefore of key importance for the struc- 
ture of the insulating phases and the phase transitions 
of the IHM. The investigation of the optical conductivity 
in the critical region is therefore a demanding task for 
future work on the complex physics of the ionic Hubbard 
model. 
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FIG. 9. Scaling behavior of the staggered bond-density 
component in the IHM near the center of the chain for [7 = 
(BI, diamonds), U — 2.1t (CI close above the transition, cir- 
cles), U = 4.0t (CI far above the transition, up triangles) and 
in the Hubbard model for U — 2.1t (squares) and U = 4. Of 
(down triangles). 
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